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EXOSAT observations of four of the five white dwarfs which 
we had obtained observing tine for were obtained in December 1985 
and February 1986. The stars observed were WD1031-115, 
WD0004+330 , WD0109-264, and WD1615-154. Counting rates were 
unexpectedly low, indicating that these objects have a 
substantial amount of x-ray absorbing matter in their 
photospheres. Analysis of these data is continuing. EXOSAT was 
not able to obtain data on the unusual magnetic white dwarf GD 
356 before it ceased operations in March 1986. 

The research on this grant took an unexpected turn with the 
discovery of 9.25 minute pulsations in the DA white dwarf V471 
Tauri. X Worked closely with the Goddard group <K. Jensen, R. 
Petre, J. Swank) on the preliminary interpretation of this 
object. A paper describing ou&" results has been accepted by the 
Astrophvsical Journal Letters and a preprint is appended to this 
report. In February 1986 a second observation of this star 
confirmed the existence of the 9.25 minute pulsations. Our 
analysis of this star is hampered by the sparseness of the grid 
of model atmospheres which is available, and a finer model grid 
will be calculated shortly. 

Another unexpected turn of events came with the availability 
of EXOSAT data on the hot white dwarf H1504+65. This unique, 
freakish object is probably the only object we know of which has 
neither hydrogen nor helium in its photosphere. I was second 
author and principal data analyst on a long paper which has been 
accepted by the Astrophysical Journal. 

In a related paper, my graduate student Peter The j 11 and I 
analyzed the existing data on the white dwarf Sirius B to see 
what constraints from other data can be placed on the properties 
of this star. EXOSAT has obtained some spectra for Sirius B. 
This paper has been accepted for publication by the Publications 
of the Astronomical Society of the Pacific. 
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Abstract 


We report the optical identification of the bright soft X-ray 
source H1504+65, originally discovered by the HEAO- 1 A-2 Low Energy 
Detectors. The optical counterpart, selected on the basis of an objec- 
tive prism plate from the Burrell Schmidt, exhibits a steep, virtually 
featureless, continuum from optical wavelengths through the ultra- 


violet. EXOSAT observations of the object confirm the identification 
and provide important soft X-ray spectral information. IUE and 
Voyager 2 spectra, combined with ground based spectroscopy, reveal the 
continuum extends without inflection out to the Lyman limit, where 
interstellar absorption cuts it off. The energy distribution is 
characterized by a single power law of the form F a from 912 & to 

A 

8000 &. At high resolution the spectrum is remarkably featureless 
except for a broad 0 VI absorption feature at X3494.7 and C IV in ab- 
sorption at 3934 R and emission at 4658 %. The presence of 0 VI is 
supported by an 0 VI resonance feature seen in the Voyager data. Weak 
Lyman alpha absorption is detected in small aperture IUE spectra, 
..=^-which is consistent- with interstellar- ab so rp t iorilP^’ ~ ~ * ------ 

The soft X-ray data, which are not consistent with any published 
models for degenerate stars having predominantly H or He compositions. 


suggest that H1504+65 is a metal rich, near degenerate object, of very 
high temperature. Qualitative considerations based on the optical 
spectra, the observed ultraviolet spectrum, and interpretations of the 
X-ray fluxes indicate that it is definitely hotter than PG1159-035 and 
similar objects. Our best estimate of the temperature is 160,000 K 
with considerable uncertainty (roughly 30,000 K) since no model or 


black body is consistent with the available data. Evolutionarily, it 
is similar to the central stars of planetary nebulae, of recent origin, 
and currently in a state of rapid thermal evolution. 
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Abstract 

We report the optical identification of the bright soft X-ray 
source H1504+65, originally discovered by the HEAO- 1 A-2 Low Energy 
Detectors. The optical counterpart, selected on the basis of an objec- 
tive prism plate from the Burrell Schmidt, exhibits a steep, virtually 
featureless, continuum from optical wavelengths through the ultra- 
violet. EXOSAT observations of the object confirm the identification 
and provide important soft X-ray spectral information. IUE and 
Voyager 2 spectra, combined with ground based spectroscopy, reveal the 
continuum extends without inflection out to the Lyman limit, where 
interstellar absorption cuts it off. The energy distribution is 

characterized by a single power law of the form F. ot A ^ from 912 & to 

A 

80Q0 &. At high resolution the spectrum is remarkably featureless 
except for a broad 0 VI absorption feature at A3494.7 and C IV in ab- 
sorption at 3934 X and emission at 4658 X. The presence of 0 VI is 
supported by an 0 VI resonance feature seen in the Voyager data. ' Weak 
Lyman alpha absorption is detected in small aperture IUE spectra, 
which is consistent with interstellar- absorption.“~ ■” 

The soft X-ray data, which are not consistent with any published 
models for degenerate stars having predominantly H or He compositions, 
suggest that H1504+65 is a metal rich, near degenerate object, of very 
high temperature. Qualitative considerations based on the optical 
spectra, the observed ultraviolet spectrum, and interpretations of the 
X-ray fluxes indicate that it is definitely hotter than PG1159-035 and 
similar objects. Our best estimate of .the temperature is 160,000 K 
with considerable uncertainty (roughly 30,000 K) since no model or 





I. INTRODUCTION 


The soft X-ray sky is a relatively recently explored portion of 
the electromagnetic spectrum ( E v 0.1-1 keV, \ v 10-100 X) . The HEAO-1 
A-2^ experiment Low Energy Detectors (LED) performed the most sensitive 

^The A-2 experiment or. HEAO-1 was a collaborative effort lead by E. 
Boldt of Goddard Space Flight Center (GSFC) and G. Garmire of Penn 
State with collaborators at GSFC, Penn State, The Jet Propulsion Labor- 
atory and the University of California, Berkeley. 

I _ ~ " _ 

all sky survey for point sources of X-rays in this energy range (Nugent 

i 

1 et al. 1983). Of the 114 sources cataloged, half were previously 

, unknown. 

I 
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We are attempting a systematic effort to identify these brightest 
sources, in the hopes that these bright objects are the closest or most 
conspicuous members of interesting classes. Unfortunately the HEAO-1 
data alone can not precisely locate the X-ray source but have an un- 
certainty in position (error box) of order 0.3 square degrees. 

H1504+65 is the seventh brightest source in the 0.25 keV band 
(0.18-.44 keV) observed by HEAO-1, but was undetected at 1 keV (0.44- 
2.8 keV) (Nugent et al. 1983). Objective prism plates taken with the 
0.6 m Burrell Schmidt at Kitt Peak of the H1504+65 error box revealed 
a faint (m^ v 16), extremely blue continuum object at (a, 5 = (1950.0) 

« (15^ 01 m 23.1 s , +66® 24 T 09"). Moderate resolution spectroscopy 
using the Palomar 1.5 m reflector and the Penn State 1.6 m reflector 
confirmed the blue continuum, but found no evidence for emission or 
absorption lines. High resolution optical spectroscopy using the 2.1m 


Steward Observatory reflector have placed stringent limits on the 
presence of absorption or emission features from 4000 to 5000 A, but 
did show some features from C IV and 0 VI. Ultraviolet spectroscopy 
from the International Ultraviolet Explorer CIUE) also found an ex- 
tremely blue continuum, and small-aperture spectroscopy indicated the 
presence of interstellar Lyman alpha. Voyager observations have shown 
that the continuous slope extends without break, to the Lyman limit, and 
also demonstrated the presence of 0 VI 1032. EXOSAT observations, 
using the Channel Multiplier Array CCMA) , have confirmed the identifi- 
cation of the X-ray source with the optical candidate and provided 
important soft X-ray spectral information by using several filters. 

In the following sections we describe the optical observations 
(§11), the ultraviolet observations (5III), the X-ray observations 
(§IV), interpret these observations (SV) and discuss our understanding 
of the nature of H15Q4+65 (IVIl. 


II, OPTICAL OBSERVATIONS 


a) Objective Prism 

An objective prism Illa-J plate exposure of 75° was taken of the 
H1504+65 error box on 1982 April 24, using the 0.6 m Burrell Schmidt 
with thin prism at Kitt Peak National Observatory by A. Wasilewski as 
part of a Visiting Astronomer program with J.A.N. and D. Weedman. Care- 
ful search of the H1504+65 error box by Dr. Weedman revealed a dramati- 
cally blue continuum object at a, 5(1950) ■ 15^01 m 23fl, +66°24 , 09'’ 

(Figure 1) . No emission or absorption features were apparent on the 
spectrum covered on the plate 0^3600-5400 8) down to a limit of approx- 
imately 150 8 equivalent width. The stellar image shows no proper mo- 
tion, with an upper limit of 0.06 M yr ^ (Gidas, private communication). 

b) Low Resolution Spectroscopy , 

Optical spectroscopy of the blue candidate was performed at the 1.6 
m reflector of the Black Moshannon Observatory (BMO) of Penn State by 
JAN and L. Takalo and at the 1.6m reflector of the Palomar Observatory 
by SHP and F. Marshall. The BMO observations on 1983 July 13 using a 
fiber coupled SIT spectrograph, (Ramsey, Nations and B ard^l98 1) total- 
led 3000 seconds and covered the spectral range 4500-8100 8 with a reso- 
lution of 25 8. The Palomar observations on 1982 May 25 employed a SIT 
Cassegrain spectrograph, totalled 1800 seconds and covered the spectral 
range 3500-6500 8 at a resolution of 15 8. 

These observations found a steep, blue, featureless continuum 
(Figure 2) with a magnitude of 16.1 at 5000 8. In particular the upper 
limit to the equivalent width of any absorption or emission is 6 8 at Ha and 
5 8 at Hg. 

c) High Resolution Spectroscopy- 

Higher resolution observations of H1504+65 have been obtained by 


J.B.H. and J.L. with the Steward Observatory 2.3 m reflector and the 
Multiple Mirror Telescope. Initial 2.3 ta observations were performed 
on 1985 February 12 using the photon counting Reticon detector system 
with an ultraviolet - sensitive image tube. An 832 line mm” 1 grating 
was employed which achieved an effective resolution of ^2.25 X over 
the 3900 X to 5000 X region. A total effective observation time of 

Q 

1320 s yielded a steeply rising continuum with no apparent features. 

Q 

It Is worth emphasizing that, over the optical spectral range of 
4000-5000 X, this star more closely resembles a continuous spectrum 
than any other star we have observed at this spectral resolution with 
data of similar quality, including at least 1000 PG Survey stars. Only 
DC white dwarfs cooler than 12,000 K exhibit fewer apparent features. 

Of course some central stars of planetary nebulae are assigned "con- 
tinuous" spectra but generally the quality of the observations is not 
as good with the stellar spectrum being complicated by overlying 
nebular emission. 



In particular there was no evidence for either Hy or H6, the upper 
limit for the former being 0.3 X. 


A second high resolution spectrum was obtained with the Arizona/ 
Smithsonian 4.5 m Multiple Mirror Telescope on 31 March CU.T.). The 
use of a similar 832 line mm grating afforded coverage of 4050- 
4950 X at 1 X resolution. The detector system consisted of a blue- 
sensitive image tube and jphoton counting Reticon. However, the larger 

telescope aperture and better sensitivity of the detector near 5000 X 

1 

allowed better coverage of the 4500-5000 X region in particular than 
for the 2.3 m scan. In agreement with the former, there was no con- 


elusive evidence for broad or narrow stellar absorption features. 

Upper limits for Hy, HS, He II X4686 and He II X4540 are 0.5 X or less 
(Fig. 3a). In particular, the "absorption trough" due to blends of 
C IV X4658 and He II X4686, so characteristic of the so-called PG1159 

/■t 

spectroscopic group (Weseaael, Green and Liebert 1985) is absent. 

9 

The only conceivable absorption hare consists of extremely broad, 
shallow line wings with virtually no cores. Such a possible feature 
appears near 4500 X as well, but this has no logical identification. 

» 

convincingly in both independently-redujced array spectra. This line 
appears marginally resolved, with a FWKM width of about 2.8 X. The 
corresponding C IV X4441 emission line also seen in some PG1159 stars 

I 

i 

is at best marginally present in this spectrum. 

A follow up 2.3 m observation of H1504+65 was performed on 28 May 
1985 with the same instrumentation but with the grating tilted so as to 
cover the 3200 to 4300 X region. This spectrum, with an effective 
integration time of 1380 s, is shown in Figure 3b. A broad absorption 
feature exhibiting a narrow 'central emission reversal is clearly - ' 
apparent. We identify this feature as due to 0 VI X3434.7 and note the 
possible presence of weaker features due to the same ion at 3314 X and 
3622 X. A broad but weak C IV feature is also evident at 3934 X as 
well as unidentified absorption blends blueward of that position. 

Using the CIV 4658 X emission feature as a fiducial we determine 
a heliocentric radial velocity of 0+15 km/s for H1504+65. 


III. ULTRAVIOLET OBSERVATIONS 


a) IUE 

Despite its faint optical flux CV ^ 16.24) the presence of X-ray 
flux and the extremely blue optical continuum encouraged us to request 
ultraviolet observations of the candidate using IUE CBoggess et al. 

1978). Low dispersion spectra Cresolution v 6 R) were obtained by 
J.A.N. on 1983 June 13 and 1984 October 31 with both the long and short 
wavelength cameras, covering the wavelength range from 1150-3200 R. 

The ultraviolet spectra show no sign of variability between obser- 
vations. As in the optical, the ultraviolet continuum is steeply, and 
smoothly rising toward short wavelength CFigure 4) . Because of the 
broad geocoronal emission, the standard, large aperture spectra provide 
only rather weak constraints on the presence or absence of Lyman alpha. 
Consequently i a small aperture spectrum, SWP 25323, was obtained by 
H.L.S. on February 26, 1985, revealing a modest Lyman alpha absorption 
feature. Because the sensitivity of the IUE spectrograph falls steeply 
shortward of Lyman alpha, uncertainties in continuum placement limit 

.our .ability to measure the equivalent width of this -feature . __0ur_ bes t 

estimate is that the equivalent width is 5.7 + 1 A. 

Except at Lyman alpha there are no emission or absorption features 
in the ultraviolet exceeding 20% of the continuum. Possible weak ab- 
sorption features may be present near 1258, 1290, and 1300 R with 
equivalent widths near 2 R. No sign of interstellar extinction is 
found near 2200 R. 

b) Voyager 2 

In addition to IUE , far UV observations of H1504+65 were also ob- 
tained with the ultraviolet spectrometers on board Voyager 2 


(Broadfoot et al . 1977 and 1981) . These observations covered the 
period 1984 November 19 to 1984 December 28. During this period 7684 
s of high rate data (3*84 s/spectra) and 9215 s of low rate (575 s/ 
spectra) data were obtained with the source near the center of the 
field of view. The data were reduced in a manner similar to that 
described in Drilling, Holberg, and Schoenbemer (1984) and employed 
the calibration of Holberg jit al . ( 1982 ) . 

Resultant 912 to 1083 A fluxes are shown in Figure 2. These data, 
which are obtained at 9.26 A per channel, have an effective spectral 
resolution of ^ 25 £. At t* s relatively low tfssolution a statistical- 
ly significant absorption feature centered on a wavelength of Il032 A 

i 

is apparent. In light of the presence of 0 VI in the visual spectra we 
feel secure in identifying this feature as a blend of the resonant . 

0 VI doublet lines XX 1031.9 and 1037.6. A weaker, but still signifi- 
cant feature, is also present at v 977 £. 


IV. 


X-RAY OBSERVATIONS 


a) HEAO-1 A-2 Low Energy Detectors (LED's) 

H1504+65 was discovered by the HEAO-1 LED experiment as a point 

source of X-rays in the 0. 18-0.44 keV band. The HEAO-1 LED’s were 

collimated thin window proportional counters which scanned along lines 

of constant ecliptic latitude once every 33 minutes, and stepped along 

the ecliptic equator at a rate of 1° per day (Rothschild et al. 1979). 

In this way virtually the entire sky was surveyed for X-ray sources 

with a positional resolution of approximately 0.3° x 1°. 

HEAO-l scanned, H1504+65 in 1977 December and 1978 May. The 

0.18-0.44 keV X-ray intensity in the two observations showed no sign 

of variability and were consistent with a count rate of 5.9 + 1.4 

counts s 1 . No source was detected at this position in the 0.44-2.8 

keV band down to a 2a limit of 1.7 counts s 

The spectrum of pulse heights induced in the LEDs by H1504+65 was 

very soft, confined to the lowest three pulse height channels. Within 

photon counting statistics the spectral shape is the same as that 

observed in HZ 43, a very low energy black body spectrum, T < 10^ K, 

with no indication of absorbing gas. 

The determination of absolute fluxes from the HEAO-l data is 

rather model dependent. For illustrative purposes, we have determined 

the flux in the 0.18-0.44 keV band for a 150,000 K black body input 

20 -2 

spectrum and a hydrogen column density of 10 cm , to be 7.4 + 1.8 x 
-12 2 -1 

10 erg cm s . However we note that since the spectrum is very 
soft, this number should be regarded with care in any future atmospher- 
ic modeling efforts. .The need for caution in modeling the X-ray 
counting rates for sources? of this type has been known for some time 


(Broadfoot et al. 1977 and 1981) . These observations covered the 
period 1984 November 19 to 1984 December 28. During this period 7684 
s of high rate data (3.84 s/spectra) and 9215 s of low rate (575 s/ 
spectra) data were obtained with the source near the center of the 
field of view. The data were reduced in a manner similar to that 
described in Drilling, Holberg, and Schoenbemer (1984) and employed 
the calibration of Holberg e£ al . (1982) . 

Resultant 912 to 1083 A fluxes are shown in Figure 2. These data, 
which are obtained at 9.26 A per channel, have an effective spectral 
resolution of v 25 8. At this relatively low resolution a statistical- 
ly significant absorption feature centered on a wavelength of 1032 A 
is apparent. In light of the presence of 0 VI in the visual spectra we 
feel secure in identifying this feature as a blend of the resonant , 

0 VI doublet lines XX 1031.9 and 1037.6. A weaker, but still signifi- 
cant feature, is also present at v 977 


(Shipman 1976, Kahn et al. 1984). 

b) EXOSAT Low Energy Experiment CMA 

The EXOSAT Low Energy Experiment (LE) observed H1504+65 on 1984 
March 28 using a channel multiplier array (CMA) at the focus of an 
imaging telescope (de Korte et al. 1981) . It detected a bright X-ray 
source coincident, within the 10” error radius, with the blue optical 
candidate discussed in §11, confirming the X-ray source identification. 
The CMA has no intrinsic spectroscopic capability, but observations 
were conducted through three different filters, allowing inference of 
spectral shape by the variation in observed intensity. The filters 
used were the thin Lexan (Lexan 3000) , Parylene-N with aluminum (Al) , 
and boron on polypropylene (B). The observed rates from H1504+65 were 
6.97 + .04, 0.47 + .02 and 0.017 + .002 counts s \ for the Lexan 3000, 
Al and B filters respectively. 

FFT analysis of the thin Lexan photon arrival times showed no 
evidence for periodicity. The upper limit to variation is 5 % for 
timescales between one second and one hour. Rates seen in the other 
bands similarly show no evidence of periodicity. 


V. INTERPRETATION 


While there is a considerable amount of data available on this 
star, there is no single model in the literature which fits it all, 
and the discrepancies between model and observation, particularly in 
the X-ray spectral region, are extreme rather than marginal. This 
suggests that the chemical composition of H1504+65 is rather unlike the 
composition of existing models Ci.e., pure H, pure He, or solar 
composition), though further modeling would be necessary in order to 
demonstrate this suggestion conclusively. Our inability to fit models 
to data makes it necessary to appeal occasionally to qualitative or 
semi-quantitative arguments. In this section, we discuss the inter- 
pretation of the data, wavelength region, and in section VI we discuss 
our overall interpretation of this highly unusual object. 

a) Optical Observations 

The most striking feature of the optical spectra is the absence of 
H I, He II lines, lines which are. found in all other degenerates and 
subdwarfs, and in planetary nebula nuclei with sufficiently high 
quality spectr a. T he upp er limits of 0.3 % ( for Hy andHS) and 0.5 
(for HfS, He II 4686, and He II 4540) dictate that if this object has a 
H- or He-rich atmosphere, it is considerably hotter than the highest 
temperature models in the grids of Wesemael et al. C1980) or Wesemael 
(1981) m- 200,000 (log g = 7) for the H-rich models and 150,000 for 
the He-rich models (this is the highest temperature model for which 
Wasemael quotes He line strengths). 

A purely qualitative analysis, in comparison with the PG 1159 
objects and with the hottest main sequence stars, supports the idea of 
an extreme temperature for this object. The hottest 0-type stars like 


(Shipman 1976, Kahn et al. 1984). 

b) EXOSAT Low Energy Experiment CMA 

The EXOSAT Low Energy Experiment (LE) observed H1504+65 on 1984 
March 28 using a channel multiplier array (CMA) at the focus of an 
imaging telescope (de Korte et al. 1981) . It detected a bright .X-ray 
source coincident, within the 10" error radius, with the blue optical 
candidate discussed in §11, confirming the X-ray source identification. 
The CMA has no intrinsic spectroscopic capability, but observations 
were conducted through three different filters, allowing inference of 
spectral shape by the variation in observed intensity. The filters 
used were the thin Lexan (Lexan 3000) , Parylene-N with aluminum (Al) , 
and boron on polypropylene (B) . The observed rates from H1504+65 were 
6.97 + .04, 0.47 + .02 and 0.017 + .002 counts s 1 , for the Lexan 3000, 
Al and B filters respectively. 

FFT analysis of the thin Lexan photon arrival times showed no 
evidence for periodicity. The upper limit to variation is 5% for 
timescales between one second and one hour. Rates seen, in the other 
bands similarly show no evidence of periodic ity . - . ^ 


Zeta Puppis, with TCeff) ■ 46,000 + 4,000 K (de Jager 1980) have both 
C III and C IV lines, as do the hottest Wolf-Rayet stars which show 
carbon. H1504 has only C IV lines. This type of ionization argument 
is relatively independent of chemical composition, and hence does not 
require that we assume a R- or He-dominated atmosphere. 

A closer, and more instructive, comparison can be made between 
H1504 and the PG 1159 objects. Sion, Liebert, and Starrfield (1985) 
have discussed similar 0 VI and C IV features in five members of this 
class of objects, including PG 1159-035 itself. In the Visual the 
spectra of PG 1159 objects are characterized by steep, nearly Rayleigh- 
Jeans continua with broad features due to He II, C IV, possibly Si IV, 
and no H I Balmer lines (Wesemael, Green, and Liebert 1985) . In the 
optical ultraviolet Si.on, Liebert, and Starrfield find that 0 VI 3434.7 
and C IV 3934 are strong features in the PG 1159 objects they examined. 
Both are evident in our data, but they are substantially weaker in 
H1504. Further, the C IV 3689 line, strong in the PG 1159 objects, is 
absent in H1504. Because it comes from a different multiplet than the 
C IV 3934 line, the difference in the 3934/3689 line ratios (> 1 in 
PG 1159 objects and < 1 in H1504) is not a definitive indication of a 
misidentification of the C IV 3934 line, but it is an indication of the 
need for further modeling and investigation of the C spectrum. The 
weakness of the C IV and 0 VI lines which we do detect could conceiv- 
ably be due to an underabundance of 0 VI and C IV relative to the PG 
objects, but is is more likely due to a significantly higher tempera- 
ture. Additionally, the presence of He II in PG 1159 and its absence 
in H15Q4 suggests that if the atmosphere of H1504 is He dominated, or 
has a similar He content to that of PG 1159, that H1504 is considerably 


hotter so that He is more ionized (to He III) . 

The similarity between K1504 and the PG objects does not extend to 
the IUE spectral region, but again a higher temperature for H1504 
seems to be the cause. Other dissimilarities (the absence of C IV 1550 
and N V 1240) could be either temperature or abundance effects. 

The 0 VI 3434.7 line found in H1504 is at least as strong as found 
in PG 1159 objects. Note that the width of this profile is also at 
least as great as seen in PG 1159.. While here as with other data a 
model atmosphere analysis will be required for definitive results, this 
profile similarity suggests some statements about surface gravity. 
Wesemael, Green and Liebert C1985) estimated log g - 7 for the PG 1159 
group based on comparison of the He II lines with models. We believe ' 
that the similarity in the 0 VI profile of H1504+65 indicates a cor- 
respondingly high surface gravity. 

b) Ultraviolet Continuum 

Figure 2 demonstrates that the continuum rises very steeply into 
the ultraviolet. The SWP data alone, as well as the SWP data plus the 

visual point, fit a relationship of H^ ot % 4. 1+0. 1^. w ^ich the same 

2 i " ' 

as H a v *' ). The LWP and Voyager data are both noisier and are some- 
what offset from the SWP data; in our view, they are consistent with a 
power law with a spectral index of 4.1 + 0.1. (Since IUE fluxes are 
generally quoted in wavelength units, we define a spectral index x by 

“X 

the equation H^ a % . 

One can take two approaches to interpreting these continuum 
slopes. Calibration uncertainties are the smallest if one only con- 
siders the data from a particular spectrograph on a particular instru- 
ment, and the one with the greatest leverage and quality here is the 


Zeta Puppis, with TCeff) * 46,000 + 4,000 K C<ie Jager 1980) have both 
C III and C IV lines, as do the hottest Wolf-Rayet stars which show 
carbon. H1504 has only C IV lines. This type of ionization argument 
is relatively independent of chemical composition, and hence does not 
require that we assume a H- or Ke-dominated atmosphere. 

A closer, and more instructive, comparison can be made between 
H1504 and the PG 1159 objects. Sion, Liebert, and Starrfield (1985) 
have discussed similar 0 VI and C IV features in five members of this 
class of objects, including PG 1159-035 itself. In the visual the 
spectra of PG 1159 objects are characterized by steep, nearly Rayleigh- 
Jeans continua with broad features due to He II, C IV, possibly Si IV, 
and no H I Balmer lines (Wesemael, Green, and Liebert 1985). In the 
optical ultraviolet Si.on, Liebert, and Starrfield find that 0 VI 3434.7 
and C IV 3934 are strong features in the PG 1159 objects they examined. 
Both are evident in our data, but they are substantially weaker in 
H1504, Further, the C IV 3689 line, strong in the PG 1159 objects, is 
absent in H1504. Because it comes from a different multiplet than the 
- C IV 3934 line, the difference in .the 3934/3689 line ratios . (>. l_in 
PG 1159 objects and < 1 in H1504) is not a definitive indication of a 
misidentification of the C IV 3934 line, but it is an indication of the 
need for further modeling at,& investigation of the C spectrum. The 
weakness of the C IV and 0 VI lines which we do detect could conceiv- 
ably be due to an underabundance of 0 VI and C IV relative to the PG 
objects, but is is more likely due to a significantly higher tempera- 
ture. Additionally, the presence of He II in PG 1159 and its absence 
in H1504 suggests that if the atmosphere of H1504 is He dominated, or 
has a similar He content to that of PG 1159, that H15Q4 is considerably 


IUE SWP spectrograph.. Fig. 6 compares our observed value to those 
given by models; the uncertainty in our observed value includes equal 
contributions of 5% from photon statistics and 5 % from the uncertainty 
in the IUE calibration. It is clear that the observed "super- infinite" 
spectral index is, even considering the substantial uncertainties from 
calibration, higher than those produced by the H-rich models of 
Wesemael et al. (1980), and is hotter chan the Wesemael l s (1981) 
hottest He model with T(eff) * 150,000 K. A blackbody "fit" gives a 
nominal 2- sigma lower limit to the temperature of 500,000 K; this limit 
should not be taken too seriously. 

Anojther approach is to use the widest baseline possible, where 

I 

good data exists. Comparison of the Vovager and optical data gives a 
spectral; index of 4.3; comparison of the IUE SWP data and the optical 
gives a spectral index of 4.0. UV/optical spectral indices from the 
Wesemael models are about 0.08 higher than the UV indices plotted in 
Fig. 6; one comes up with a similar conclusion, that neither the H- 
rich nor the He-rich Wesemael models fit the data, that the star is 
very hot, and that it is hotter than the PG 1159 objects. 

e) Lyman Alpha 

The small-aperture SWP spectrum did show a Lyman alpha feature. 
This is clearly interstellar. If it were stellar, the indicated 

3 

effective temperature from the Wesemael models would be (80-100) x 10 
K, not an implausible value, but then the equivalent width of HB or Hy 
would be at least 1.8 A for a continuum 30 A from line center. A 
relationship between the equivalent width of interstellar Lyman alpha 
and the interstellar column density N is given by Morton (1967) as 

Cl 


*« / i \ 


■* A m 


,.- 10 r „ . 1/2 


While the rapidly changing sensitivity of the IUE , and the somewhat 

noisy small aperture spectrum, make measurement of W uncertain, the 

observed value of W - 5.7 + 1 A leads to an interstellar column 

20 —2 

density of (.0.6 + 0.2) x 10 atoms cm , While the fits to the X-ray 
data reveal no clear choice of model or Njj Csee below) , this value is 
reasonably consistent with the X-ray data, particularly the very soft 
data point from EXOSAT . 

d) HEAO-1 

There are a fair variety of approaches to analyzing the X-ray 
emission from individual white dwarf stars (see, e.g. , Shipman 1976, 
Auer and Shipman 19.77, Malina, Basri, and Bowyer 1982, Martin et al. 
1982, Kahn et al. 1984, Holberg, Wesemael, and Huben? 1984, Petre, 
Shipman, and Canizares 1985). The difficulty faced by all these 
authors is that the number of spectrally independent data points in 
the X-ray band is typically one or, at most a few, low resolution 
measurements of the integrated flux. Balancing this difficulty is the 
tremendous leverage offered by the extremely steep dependence of the 
X-ray. flux on the physical parameters of interest , .usually _T(eff)_.and 
n(He) , and its virtual independence of surface gravity. 

We face an additional difficulty in analyzing the HEAO-1 data in 
that the detector calibration is a strong function of the shape of the 
assumed incident spectrum. For ultrasoft sources like white dwarf 
stars, the convolution of the steeply falling incident spectrum with 
the sharply rising detector response means that virtually all the 
counts come from the softest energy' which the instrument can detect. 
This means that the distribution observed in the detector is unchanged 
over a wide range of assumed incident spectra, but the detector 
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efficiency changes dramatically. For example, the soft X-ray flux 
inferred from a given count rate changes by a factor of twenty when 
the assumed incident spectrum is changed from a 100,000 K black body to 
a 200,000 K black body. Hence we can not infer a temperature or even a 
flux in the soft X-ray band without knowledge from other data of the 
incident spectral shape. 

To escape this extreme model dependency we have chosen to form a 
completely empirical quantity, IL /l, ] , based on the ratio of X-ray to 
optical flux between H1504 and HZ 43. We define 


[L /L ] 
x v J 


If x /f v ] H1504+65 

^ f x /f v^HZ43 


ft 
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where the f's are the observed fluxes. To first order the ratio f 
takes out differences in distance and stellar radius, so [L /L ] should 

X V 

be one, if H1504 and HZ 43 have the same spectrum CHZ 43 was selected 
as a reference because it was well observed optically and by HEAO-1 and 
EXOSAT ) . [L /L ] is found to be 8.1 + 2, indicating that H1504 has a 
spectrum proportionately emitting more X-rays than HZ 43. 

To achieve quantitative temperature estimates it is necessary to 


invoke models, but by making the comparisons in terms of [L^/L^] we 
hope that model inadequacies will also cancel to first order. We used 
a 60,000 K pure hydrogen model to represent the spectrum of HZ 43. 
Uncertainties in the temperature or spectrum of HZ 43 affect this 
scheme, but by a small amount; a 3000 K change in the temperature of 
HZ 43 changes the value of [L /L ] by only 1.5, a small effect in light 

X V 

of the tremendous leverage offered by observations of thermal X-rays. 

An additional correction to be made is that for interstellar 
absorption. Nugent et al. Q.983) in Fig. 4, give corrections which 


can be applied to black bodies. In those cases where we tried to fit 
models to the data, we used the black body temperature which came 
closest to fitting the model spectrum in the 100-200 eV region. 

Now consider fitting the data for various possible chemical com- 
positions. If H1504+65 were H-rich, like HZ 43, the spectrum in the 
100-200 eV region would resemble that of a 300,000 black body, for the 
temperatures in the vicinity of 10^ K which we derive. Correction for 

interstellar absorption indicates that fL /L ] , its value at the 
r x v o 

surface of H1504, is 11+4. (This uncertainty includes the uncertain- 
ty in the zero point caused by uncertainties in what parameters to use 
for HZ 43.) This would fit a pure H model at TCeff) = 98,000 K. 

Adding very substantial quantities of trace He, with N(He) « 0.01, 
characteristic of the DO stars, would increase the derived temperature 
to 150,000 K. 

The He-rich models of Wesemael (1981), near T(eff) - 150,000 K, 

resemble 400,000 K black bodies, and so the correction for interstellar 

absorption is a bit less than in the K-rich case, with fL /L ] = 

.X V o 

-9. 6 + 4. .There is some modest gravity sensitivity .to the value -.of^— 

fL /LI determined from the models, and so the derived temperature 
X v 

using He-rich models is 156,000 K for log g = 8, with an estimated 
uncertainty of 10,000 K produced by an uncertainty in log g. The grid 
of models is rather sparse; with IL^/L^] increasing by a factor of 2 
when one decreases log g from 8 to 7 at 200,000 K, and with a change 
in the other direction at 100,000 BC. 

The black body analysis proceeds along similar lines, where one 
used Nugent et al.’s Fig. 4 to estimate corrections for interstellar 
H (which vary with the blackbody temperature) , and one compares a 
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the assumed incident spectrum is changed from a 100,000 K black body to 
a 200,000 K black body. Hence we can not infer a temperature or even a 
flux in the soft X-ray band without knowledge from other data of the 
incident spectral shape. 

To escape this extreme model dependency we haye chosen to form a 
completely empirical quantity, IL /L ], based on the ratio of X-ray to 
optical flux between H1504 and HZ 43. We define 


[L /L 1 


If x /f v 1 H1504+65 


x v J [f /f 

1 x v HZ43 


where the f’s are the observed fluxes. To first order the ratio f /f 

x' v 

takes out differences in distance and stellar radius, so [L /L ] should 

> x v J 

be one, if H15Q4 and HZ 43 have the same spectrum (HZ 43 was selected 
as a reference because it was well observed optically and by HEAO-1 and 
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spectrum proportionately emitting more X-rays than HZ 43. 

To achieve quantitative temperature estimates it is necessary to 


invoke models, but by making the comparisons in terms of [L^/L ] we 
hope that model inadequacies will also cancel to first order. We used 
a 60,000 K pure hydrogen model to represent the spectrum of HZ 43. 
Uncertainties in the temperature or spectrum of HZ 43 affect this 
scheme, but by a small amount; a 3000 K change in the temperature of 
HZ 43 changes the value of by only 1.5, a small effect in light 

of the tremendous leverage offered by observations of thermal X-rays. 

An additional correction to be made is that for interstellar 
absorption. Nugent et al. (19.83) in Fig. 4, give corrections which 


black body spectrum to the HZ 43 spectrum and calculates IL /L ] . 

X VO 

Since the black body spectrum is considerably softer than the spectra 
of H or He models at similar temperatures, the correction for inter- 
stellar absorption absorption is greater. The best fit is for T * 

(166 + 10) x 10 3 K. 

e) Attempts to Model the EXOSAT Count Rates 

The EXOSAT observations of H1504+65 consist of count rates ob- 
tained from three overlapping, but significantly different broad band 
filters covering the wavelength range 100-240 R. At the very lowest 
energies the response in the EXOSAT filter bands is a more senitive 
measure than the HEAO-1 proportional counter. Potentially then, we can 
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constrain the effective temperature of H1504+65 in terms of published 
models. The approach is as follows. As we remain largely ignorant of 
the true atmosphere composition of H1504+65 we investigate four quite 
different types of model atmospheres; a simple black body, pure H 
models (Wesemael _et al. 1980) , pure He models (Wesemael 1981) , and 
solar composition models characteristic of the central stars of 
planetary nebulae (Hummer and Mihalas 1970). For each model (except 
black body), we employ gravities in the range log g = 7.5 to 8.0 and 
normalize to the observed visual magnitude (V = 16.24). EXOSAT count 
rates are then computed as a function of temperature and interstellar 
column density, N . For each filter the observed count rate and its 
associated statistical uncertainly correspond to a band of temperatures 

= Ng plane. Acceptable models correspond to 
regions where the xhree filter defined bands mutually overlap. 

No one model, or set of models, simultaneously predicts the 
observed EXOSAT count rates in all three EXOSAT filters. More over. 
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very few models succeed in predicting observed count rates for 
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reasonable values of-N„ (< 10 ’ cm' ), In general the disagreements 

CL 

are very large. Expressed in terms of the la error bars in the count 
rates given in sec. IVb the discrepancies are at the 10-20 a level. 
Inconsistencies of this magnitude greatly exceed those which might be 
attributed to our knowledge of the effective areas used to define the 
filter band passes. Clearly we have not correctly modeled the soft 
X-ray flux from H1504+65. This is perhaps not unexpected since none 
of the models considered should be expected to provide an adequate 
representation of soft X-ray fluxes from H1504+65. A black body is of 
course a poor approximation in the soft X-ray (Shipman 1976) . The pure 
H models completely lack important He II and metal opacities and are 
consequently too bright. He rich models, which might apply to PG 1159 
objects, too thoroughly quench the soft X-ray flux, and produce almost 
no flux below the He II ionization edge at 227 % (which is where 
virtually all the EXOSAT counts are observed.) At exceedingly high 
temperatures (T >> 150,000) He rich models may begin to provide 
~suf f icientrf lux7~but~ the predicted ~r atib of ^EXOS AT ~~f i Iter rates'Ts 
inconsistent with the data. At intermediate He content the unblanketed, 
LTE, solar composition models of Hummer and 'Mihalas also produce 
mutually inconsistent EXOSAT count rates. This particular situation is 
illustrated in Figure 7 where effective temperatures and interstellar 
column densities corresponding to the observed count rates in the 
three EXOSAT filters are shown. 

It is clear that until representative model atmospheres become 
available soft X— ray fluxes from such objects will remain difficult to 
interpret. To illustrate, the EXOSAT - count rates for H1504+65 are 
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observed EXOSAT count rates in all three EXOSAT filters. More over, 



contrasted in Table 1 with those of two recognized PG 1159 objects, 

PG 1159-035 and Kl-16, and with the DA white dwarf HZ 43. The last 
column of Table 1 gives a soft X-ray to optical luminosity ratio, 
normalized to 1.0 for PG 1159-035, for all four objects. Note that 
PG 1159-035 and Kl-16 have similar luminosity ratios while H1504+65 is 
considerably larger. 


VI. DISCUSSION 


a) What is H1504+65? 

There is considerable data for this star, but the above sections 
indicate that it does not seem easy to fit it with existing models. We 
can conclude that it is definitely not a DA white dwarf star, a hot 
analog of HZ 43. The implied X-ray temperatures are far too cool to be 
consistent with the absence of Balmer lines. Even a model with sub- 
stantial trace He, with TCeff) = 150,000 K, produces an Hy line with a 
strength, depending on gravity, of 1.2 - 2 A, incompatible with the 
optical spectrum. The disagreement between model and data is 
sufficiently great that we believe we can rule out this explanation. 

Our attempts to fit the Wesemael pure He models to this object 
also did not succeed, but the discrepancies between model and observa- 
tion are not as drastic. The X-ray data suggest a temperature of 
160,000 K. Wesemael does not provide any calculations of He line 
strengths for temperatures higher than 150,000 K Csee his Figure 7). 

The models indicate that the He lines are weakening with increasing 
^temperature,- but that they are still quite strong at 150,000 K..(l. 26. A 
for log g ■ 6 and 3.23 A for log g - 8). Additional modeling would be 
required to decide, definitively, whether the He lines fall off fast 
enough' with increasing TCeff} so that a Ke-dominated atmosphere with 
a TCeff} of around 160,000 K would be consistent with an upper limit 
to the strength of this line of 0.5 A, Additional problems with 
the He interpretation are the steepness of the ultraviolet- continuum 
and the failure to fit the various EXOSAT data points. 

If the primary constituent of the atmosphere is neither H nor He, 
what could it be? The only spectral features that we see are C IV and 
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0 VI. While any definitive investigation of the possibility that we 
are seeing a carbon-oxygen dominated atmosphere must await further 
modeling, some qualitative apd semi-quantitative considerations indi- 
cate that this suggestion is reasonably consistent with the features 
that we observe. At this high tempeTautre, with an electron density of 
10^ (characteristic of optical depth unity), C is mostly C V, with 
some C IV at shallower depths. C V, a He— like ion, is essentially 
unobservable; its ionization potential is very high, 391 eV, and the 
resonance lines are in the X-ray part of the spectrum. The only 
possible problem with the C spectrum that we observe is whether the 
weak C IV lines that are seen in the optical spectrum are consistent 
with the IUE data, which shows the C TV resonance lines to be marginal 
at best. In a high-density, high-temperature atmosphere, 0 will exist 
as 0 V, 0 VI, and 0 VIII. Except for the 0 VI resonance doublet at 
1031 and 1037 A, which we have detected, these species all have their 
resonance lines in the EUV spectral region, where we couldn't see them. 
Whether the features that we detect (0 VI 1032, 1038, and 3434) are 
consistent with absence of the 0 VI lines at 3811 and 3834 A remains to 
be seen. Perhaps the 3811 and 3834 A lines are heavily broadened, or 
even shifted to become the unidentified features seen at 3860 and 3885A. 

To summarize, the line spectrum that we see suggests a C-0 
dominated atmosphere, with no obvious resonance lines missing, although 
the absence of C IV 1548, 1550 and 0 VI 3811, 3834 is perhaps puzzling. 
Without detailed models, we are not prepared to claim that we have dis- 
covered a white dwarf star of this unique, if not unexpected, composi- 
tion. Such modeling is clearly a high priority task for the immediate 
future, as is the calculation of models with the composition and 


T ._ > 150,000 K. It Is nevertheless true that, unlike the cooler hot 
eff 

degenerates like HZ A3 and the DO stars, existing models are unable to 
account for many pieces of information: the different X-ray count 
rates from EXOSAT, the ultraviolet continuum slope, and the absence of 
II I and He II lines, although it is possible that a He-dominated model 
could be generated which would produce the necessary temperature and 
the absence of He II, 

So far, very little has been said about surface gravity, except 
for some comparisons with the PG 1159 objects which indicate that 
H1504+65 has a surface gravity comparable to that of PG 1159 
(log g ■ 7) . There are, however, independent arguments, based on its 
faintness and position in the sky, which constrain it to be a degenerate 
or near degenerate. The optical flux from a stellar photosphere is 
relatively independent of composition, varying by about 30% from H- to 
He-dominated models. The fundamental relationship between monochro- 
matic Eddington flux H^, the flux received at earth f^, the stellar 
radius R, and its distance D is CShipman 1979) 
f = 4ir H CR/D) 2 

y v _ . 

With H1504+65 at a galactic latitude of 46°, one can use plausibility 

arguments to constrain acceptable values of R. A temperature of 
160,000 K determines the value of H in eq. 1. From the observed 
value of f y , we have relationship between R and z: (R/0.03 Rq) = 

(z/850 pc). (R » 0.03 Rg corresponds to log g = 7, for a plausible 
mass of 0.6 solar masses.) 850 pc is already implausibly large as a 
distance from the galactic plane; gravities less than log g = 7 are 
probably excluded simply because of this object’s faintness and its 
position in the sky. We conclude that H1504+65 is a degenerate or 
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near-degenerate, with a temperature of (160 + 30) x 10 K, and log g 
greater than 7. Our error estimate is greater than that suggested by 
any of the individual model fits, and is a rather subjective one which 
allows for the absence of any model to fit the data. 

b) Evolutionary Considerations 

Cautiously accepting our temperature of 160,000 K, and our limit 
that 7 <log g < 8, we can place H1504+65 in the HR diagram, as shown in 
Fig. 10. The evolutionary tracks are from Kawaler, Hansen, and Winget 
(1985) . The evolutionary time scales associated with these objects are 
quite short, as was dramatically confirmed when the pulsation period 
change in PG 1159-035 was observed by Winget et al. (1985) . Depending 
on the mass of H1504, the time that it would take to cool to the 
temperature c-f PG 1159-035 is only of the order of 10^ yr. Since 
[L x /L ] of H1504 is 600 times the corresponding value for PG 1159, this 
change in X-ray brightness is about 0.2 percent per year or 2 percent 
in a decade. These strictly empirical considerations suggest that such 
a change might be observable with a sufficiently long-lived spacecraft, 
especially considering that one could make a differential measurement 
in comparison with a constant X-ray source with a similar spectrum such 
as HZ 43. 

It is also apparent from Figure 2 that H1504+65 is in the part of 
the H-R diagram occupied by the hotter planetary nebula nuclei (PNNs) 
as are the PG 1159 objects. The PNN of the nebula Kl-16 also exhibits 
pulsations similar to PG 1159-035 (Grauer and Bond 1984) . Many of the 
PNNs showing strong 0 VI and C IV emission features are believed to 
have very high temperatures and He/CNO-enriched surfaces. Since at 
least most white dwarfs are believed to be descended from planetary 


nebulae, one may ask where is the nebula around H1504+65 (and likewise 
the isolated PG 1159 objects)? One possible answer is the occurrence 
of a late helium shell flash, resulting in the ejection of the remain- 
ing residual hydrogen, the exposure of helium and/or CNO-rich material 
at the surface, and the return of the PNN to the region of the Harmon- 
Seaton sequence after the original hydrogen-rich nebula has dissipated 
(Iben et al. 1983) . Such a scenario may be developed to explain any 
highly unusual surface abundances in H1504+65 in particular. 

A third evolutionary consideration is the clear difference in 
temperature between objects like H1504+65 and PG 1159-035 entering the 
white dwarf stage and the hottest of the more numerous DA white dwarfs 
with hydrogen-rich surfaces. One decade after HZ 43 and Feige 24 were 
discovered as EUV sources, it has been difficult to establish any DA 
white dwarfs which are substantially hotter than these well-analyzed 
stars near 50,000-60,000 K. While ^ 100,000 K DA stars would be 
expected to show up as strong soft X-ray sources, none were found in 
the X-ray surveys (Nugent et al. 1983) . Fleming, Liebert and Green 
(1985) find no DA star hotter than about 80,000 K among the complete 
sample of over 300 white dwarfs (and several hundred more subdwarf 0 
and B stars) in the PG Survey. Yet six PG 1159 stars near 100,000 K 
and log g 'v 7 with helium-rich atmospheres, were found from the same 
sample, and H1504+65 appears in this sample as well. 

The hottest stars entering the white dwarf sequence with hydrogen- 
rich compositions include the DAO stars with traces of helium (such as 
HZ 34; Wesemael et al. 1985) and a few PNNs (such as Abell 7 near 
75,000 K and log g 'W and NGC 7293 near 90,000 K and log g a, 6.6; 
Mendez, Kudritzki and Simon 1983). The differences in temperature (at 
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[L /L ] of H1504 is 600 times the corresponding value for PG 1159, this 
change in X-ray brightness is about 0.2 percent per year or 2 percent 
in a decade. These strictly empirical considerations suggest that such 
a change might be observable with a sufficiently long-lived spacecraft, 
especially considering that one could make a differential measurement 
in comparison with a constant X-ray source with a similar spectrum such 
as HZ 43. 

It is also apparent from Figure 2 that H1504+65 is in the part of 
the H-R diagram occupied by the hotter planetary nebula nuclei (PNNs) , 
as are the PG 1159 objects. The PNN of the nebula Kl-16 also exhibits 
pulsations similar to PG 1159-035 (Grauer and Bond 1984) . Many of the 
PNNs showing strong 0 VI and C IV emission features are believed to 
have very high temperatures and He/CNO-enriched surfaces. Since at 
least most white dwarfs are believed to be descended from planetary 
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d) Holberg, private communication 
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FIGURE CAPTIONS 


Figure 1 — Finding chart for H1504+65 (a, 6 - 15 h 01°23?l, +66°24’ 09") . 

From the Palomar Sky Survey E Plate. 

Figure 2 — Combined spectra of H1504+65. 

Figure 3a - Multiple Mirror Telescope blue-sensitive Reticon scan of 
the 4050 to 4950 % region of H1504+65, The sharp emission line 
due to C IV X4686 and He II X4658 is indicated. Upper limits for 
Hy, HB, He II X4686 and He II X4540 correspond to equivalent 
widths of 'v 0.5 X. 

Figure 3b - A 2.3 m Steward optical ultraviolet Reticon scan of the 

3200 to 4300 % region of H1504+65. The presence of a broad, j 

i 

centrally reversed, 0 VI feature at X3434.7 and possible weaker 
features due to the same ion at 3314 and 3622 R are indicated. 

i 

Also shown is a broad shallow depression corresponding to the C IV 
feature present in other PG 1159 objects. (Sion, Liebert and 
Starrfield, 1985).. 

Figure 4 — Low resolution short wavelength ultraviolet spectrum of 

H1504+65 ._ Sum of 4 SWP exposures by IUE . Lyman alpha portion of 
spectrum is taken from small aperture spectrum showing interstellar 
absorption with equivalent width of 5.7 + 1 R, 

Figure 5 — Low resolution long wavelength ultraviolet spectrum of 
H1504+65. Sum of one LWR and two LWP exposures by IUE . 

Figure 6 — The ultraviolet spectral indices, defined over the 1250- 
1800 A range, as a function of stellar temperature. Circles ; 
Wesemael C1981) pure He models. Squares ; Wesemael pure H models. 
Crosses ; the black body relation. Observed points from PG 1159 


and H1504 are shown. 


Figure 7 — Regions of - N ff model space consistent with EXOSAT 

filter observations for Hummer and Mihalas 0-970) planetary nebula 
central star models. 

Figure 8 — The HR diagram. Error bars show the position of H1504+65 
(this paper) and the PG 1155 object.s (Holberg et al. 1985) . The 
evolutionary tracks are from Kawaler et al. Numbers in boxes are 
the ages, where the zero point is set by assigning an age of 3000 
yr at log L « 3.0, in thousands of years. The large hatched 
region is the' region occupied by the nuclei of planetary nebulae 
according to Pottasch (Planetary Nebulae , 1984, fig IX-1) , 
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Temperature, radius and rotational velocity of Sirius B. 

by Peter ThejU and Harry L. Shipman. 

Physics Department, University of Delaware, 

Newark, Delaware 19716. 

Abstract. 

Based on the theoretical work, by Ostriker and Hartwick (1968), concerning the magnetic fields of 
rotating white dwarfs, and the recent determination of the effective temperature, by Holberg et al. 
(1984), we have derived the interrelationships between radius, rotational velocity and effective 
temperature for Sirius B. 

We find that i) the radii, 0.76 to 0.81 in hundredths of solar radii, predicted by between 26000 
and 28000K and a mass of 1.0 5 M q places Sirius B somewhat above the Hamad a-Salpeter mass-radius 
relation, by about 1 or 2 tr. ii) the rotational velocity of Sirius B is predicted to lie between 0 and 
600 km/s (using 2 a limits) and that only radii which place a white dwarf above the mass-radius line 
are allowed, iii) in order to further limit the rotational velocity and the radius of Sirius B better 
determinations of T eff and magnitude are needed. If the satellite HIPPARCOS can determine the 
parallax of Sirius B with 1 milliarcsecond accuracy then the errors on and the magnitude 'mil be 
the limiting factors in the future. 

Keywords: Sirius B - rotating white dwarfs • magnetic fields - Mass/Radius relation. 

L Introduction. 

The motivation for our work is twofold. "We want to check if Sirius B falls on the mass-radius (M- 
R) relation using new data from Holberg et al. (1984). Only a few objects have so far been placed 
conclusively on that relation. We also wish to present a new method for determining rotational 
velocities for cases when the standard method, of modeling rotation ally broadened lines, doesn’t work. 

Many determinations of Sirius B’s properties have been made over the years. Greenstein, Oke and 
Shipman (1971) analyzed the H q and H^ line profiles and arrived at T gff — 32000 +/- 1000K, along 
with gravity and radius: R/R q — 0.0078 +/- 0.0002. Shipman (1976) pointed out that a pure 

Hydrogen model with T = 32000K would have the low absorption coefficient in the X-ray region 
necessary to explain the X-ray observations by Mewe et al. (1975). This value of T <ff was obtained 
using unblanketed model atmospheres. If blanketed atmospheres had been used, this value of T^ ff 
would have been lowered by 2000K. Later determinations of T eff have confirmed the general 
conclusion that Sirius B is much hotter than Sirius A, but also lowered the somewhat. Bohm- 
Vitense, Dettman and Kapranidis (1979) found T^ = 26C00 +/- 1000K by fitting theoretical models to 
IUE spectra of Sirius B. Later Martin, Basri and Lampton (1982) used HEAOl A-2 data and fitted a 
28000K model finding support for a photospheric origin of the X- rays. Using the IUE data and UV 
data from Voyager 2 Holberg, Wesemael and Hubeny find support for T eff in the range 26000 *> 28000K 
with a firm upper limit at 28000K. EXOSAT data will also provide a constraint on the T #ff , although 
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uncertainties in the EXOSAT calibration are still being investigated. 

One open question is whether, given all these determinations of T <ff) Sirius B fits on the M-R 
relation. One reason Sirius B might not fit the mass- radius relation is that rotation can increase its 
radius (Ostriker and Hartwiek, 1968). While one cam analyze the rotationally broadened line profiles 
(Pilachowski and Milkey, 1984) this method is only applicable to stellar spectra with very high 
resolution - and thus only useful for a few bright stars. Another possibility is that residual nuclear 
burning could increase the radius from its zero temperature value, as has been shown by Koester and 
Schonberner (1986) for 0.6 M q stars. A similar calculation for the 1 M q case has not been done, but 
the effects of nuclear burning would be smaller, and possibly negligible, for a higher mass. 

Ostriker and Hartwiek established the relationship between radius, magnetic field and angular 
momentum for a white dwarf of mass I.OJjM^. We can use their work to get the relationships between 
radius, rotational velocity and temperature, provided that we have an upper limit to the magnetic field 
strength and have values for orbital elements, the magnitude and the parallax. We have applied this 
method to Sirius B. 

II. Method. 

From Keplers third law one can calculate the total mass of the Sirius system and use the accurately 
known fractional mass of Sirius B to find its mass. The radius oi Sirius B is found by using the 
calculated relationship between absolute flux and effective temperature (Shipman 1979, table 2) and the 
magnitude and parallax (Gatewood and Gatewood, 1978). In using the Ostriker and Hartwiek (1968) 
results we assume that Sirius B rotates as a homogeneous sphere and that there is essentially no 
magnetic field on Sirius B. We base this last assumption on inspections of high-resolution spectra of 
the H q and H^ lines (Greenstein, Oke and Shipman, 1971, unpublished spectra). 



As we used Ostriker and Hartwicks results we found that the v =0 line fell above the Hamada- 

«q 

Salpeter mass-radius relation. We assume this is because Ostriker and Hartwiek did not sufficiently 

model the envelope of white dwarfs and that if thev had done so the v =0 line would have fallen on 

■ «q 

• 12 

the Hamada Salpeter M-R relation. Thus we put v^=0 on the M-R line for a C body. We assumed 
the angular momentum of Sirius B could be modeled by using the moment of inertia of a homogeneous 
sphere. Using a more realistic mass distribution (p(r)r 2 =const,) makes only a 20% difference in the 
equatorial rotational velocity v . 

In fig.l. we plot the results of our analysis in a M-R diagram. We have drawn lines parallel to the 
Hamad a-Salpeter relation through each point determined by velocity and mass - we are thus assuming 
that the M-R lines for non-zero velocities will look like this in a small area close to the HS line. Fig.2 
is a more detailed nomogram where the relationships between T^, radius and rotational velocities can 
be seen for a 1.05M q , homogeneous, spherical white dwarf. 
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in. Uncertainties. 

The quantities that enter our calculations are of course in error to some extent * to see how 
important these errors are we both calculated the propagation of errors and introduced 1 a shifts in the 
parallax and the magnitude to see the effects. (See fig. 1.) The sizes of the relative errors 0 ^ on the 
quantities that determine the mass of Sirius B aie: 

K«,(Mb)) 2 = l**J*)t + (3*>)) 2 + KA)) 2 + (2- A)) 2 

or 

(0.0261) 2 = (0.0246) 2 + (0.0079) 2 + (0.0030) 2 + (0.0022) 2 

where it,a,i^ and P are; parallax, angular semi-major axis of the orbit, the fractional mass of Sirius B 
and the period of the orbii. The 3% error is mainly due to the error on the parallax, a contributes 
about one tenth as much. 

The error on the radius, determined from H(T eff ), it and magnitude m, is dominated by the error on 
the magnitude: 

(VP**))’ = 0-10/5 m VlM) 2 + K,,W) ! + (<-„ ,,(H(T „)/2)f 

or 

(0.0192) 2 = (0.0166) 2 + (0.0082) 2 + (0.005) 2 

where we have assumed that the relative error on H(T eff ) is 1% - a reasonable number, as will be 
"discussed in "section IV. The error on the magnitude contributes about four times as r much as does the 
error on the parallax. We have used an absolute error on the magnitude of (0.03) 2 -i-(0.02) 2 where the 
last term takes into consideration the absolute calibration error. It is obvious that an improved 
determination of the magnitude is called for. Perhaps it is impossible, at the present epoch of the 
Sirius B orbit, and with standard techniques, to improve significantly on the value for the magnitude, 
but, as Smith and Terrile (1984) point out in their paper on the ft Pictoris disk the sort of equipment 

they used is very well suited for photometric work near bright sources, such as the Sirius system. 

How will the capabilities of the HIPPARCOS satellite change matters? If it can determine parallaxes 
to within 1 milliarcsecond then the error on Mg that it introduces will 1 be the same as that from the 
semi-major axis and furthermore the uncertainty on the magnitude will overshadow that on the 

parallax by a factor of nearly 40 in its influence on the relative error of Eg. 
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TV. Discussion of the methods. 

We determine the radius R of Sirius B by using the fundamental relation between the stellar flux f, 
the monochromatic Eddington flux at the stellar surface H(T (ff ), and the stellar distance D, from 
Shipman 1979: 

(1) R = D (f/(4*H(TJ))V 2 

Because of the sensitivity of the results we discuss to relatively small uncertainties in the values of 
various input quantities, we include a considerable discussion of the above equation and the 
uncertainties therein, and ask for the readers’ patience with what is in places a pedagogical discussion. 
Earlier discussions of errors have been given ty Shipman (1979) and, with respect to the accuracy of 
the H(T #ff ) relation, by Wesemael et al. (1980). 

An advantage of this method, and of all surface brightness methods for determining a stellar radius, 
is that you do not need to know (or guess) how much energy the star is emitting at all wavelengths. 
The disadvantage is that you must be able to determine the stellar flux H accurately. But this 
problem can be circumvented or minimized by carefully choosing the wavelength or wavelength range 
over which you fit the above equation to the data. Provided that you apply the equation at a 
wavelength which is shortwaid of the Planck maximum, the dependence of H on T <ff is relatively weak 
* linear in the black-body approximation - which minimizes the uncertainty in R produced by 
uncertainties in Tff 

In this paper we apply it in the V band, since the V flux from Sirius B is known to within 0.03 
magnitudes. In contrast, Bohm-Vitense et al. applied the equation in the IUE spectral range, which 
one might think to be superior, being closer to the Planck maximum at these temperatures. However, 

“ model calculations indicate otherwise. - increases by 33%- when -you “increase T ^ from 30,OOOK~- 

to 35,000 K, while H J363 A increases by 57% (Shipman, unpublished, blanketed, non-LTE models; the 
Wesemael et al. models are virtually identical^ as is shown in table 1). Thus, a 1,000 K uncertainty 
in T eff produces an uncertainty of 6.6% in H (and 3.3% in R) when you use this method in the V 
band, and an uncertainty of 11% in H (5.5% in R) when you use the method in the IUE spectral 
range. Thir consideration, by itself, suggests that the visual band is the one to work with. 

The ultraviolet spectral band is not as useful for application of equation (1), although it is essential if 
one is to determine T eff accurately. The correction for scattered light in the IUE spectral range is 
considerable, being smallest at the shortest wavelengths, and is uncertain, since one has to use an 
existing spectroscopic instrument rather than a specially constructed area scanning photometer. In 
addition, the uncertainty in the overall IUE absolute calibration is usually quoted at 10% (see, e.g., 
R.C. Bohlin, et al., 1980 and R.C. Bohlin, 1985). Individual readers can make their own judgments 
of the accuracy of the IUE measurements of the Sirius B UV flux, but it is at least as uncertain as 
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the overall calibration uncertainty of 10%. Our examination of the figures in Bohm-Vitense et al. 
suggests that the uncertainty in measuring f generated by the need for correcting for scattered light is 
probably half of the total correction, or, at the shortest wavelengths, about 15%, producing an overall 
error in f of (10) + (15) or 18%, which corresponds to 9% in R arising from uncertainties in the flux 
measurements alone, if you apply equation (1) in the IUE spectral range. 

An additional consideration which must be applied in an investigation of this kind is how good are 
the models; how well can you trust the H(T #ff ) relation which enters equation (1) ? This point has 
been discussed earlier by, e.g., Auer and Shipman (1977), Shipman (1979b) , and most extensively by 
Wesemael et al. (1980). The ATLAS program, used for the H(T tff ) relation used here (and in 
Shipman 1979), and the Auer-Mihalas program as modified by Wesemael, used for the Wesemael et al. 
calculations, are completely independent. The uncertainties in H g405 A for a given T eff are really quite 
small (see table 2.); at T^ — 20,000 K, comparing blanketed models with blanketed models,, the 
relative change in H 540g calulated as (H ShipmM - J is +0.0059, at T = 30,000 K it 

is -0.0032 and it is -0.0061 at T^ = 50,000 K. Using unpublished, non- LTE models calculated at 
Delaware indicates a realtive change in H g405 A calculated as (H nonLTE - ®x,xE^^LTE == 

35,000 K. This scale is quite gravity independent; going from log(g) of 8 to 8.5 at 30,000 K decreases 
the, Eddington flux, in these same relative terms, by -0.0082. The uncertainties in the H(T #ff ) relation 
are thus far smaller than that of the magnitude and can safely be estimated at about 1% or less. 

V. Results. 

On Cg. 1 one sees that the temperature range that Holberg et al. arrived at corresponds to a range 
in radii of from 0.76 to 0.81 in lOOR/R^. This is above the C 12 mass-radius relation of Hamada and 
Salpeter, but if ltr shifts are introduced in magnitude, mass and parallax then Sirius B moves onto the 
M/R relation. Not until the errors on magnitude, parallax and of course T eff are reduced, as discussed 
above, can the. issue of Sirius B’s placement on or off the mass-iadius relation be -settled. --~ *“ ~ = ““^' 

From the Holberg et al. temperatures and the Ostriker and Hartwick results we find the velocity to 
lie between 0 and 600 km/s - these are the velocity limits that a 2a circle, centered at (M,R), touches. 
Furthermore, only positive rotational velocities are allowed, so only the part of the M-R plane above 
the C * line is accessible. Any temperature and mass determinations that place Sirius B, or any other 
white dwarf, below the line can thus be ruled out. 

VI. Discussion. 

Some interesting considerations, and some suggestions for current work, emerge from a consideration 
of how the astrometric data on Sirius are likely to improve in the future. The European satellite 
HlRPARCOS can improve the parallax of Sirius significantly; the current best parallax (Gatewood and 
Gatewood, 1978) has an uncertainty of 3 milliarcseconds, and the advertised HIPPARCOS capability 
should improve this by a factor of 3. Such improvement can, in principle allow us to determine the 
radius with much greater precision than we can now. In addition, at the present time the uncertainty 
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in the mast is dominated by the uncertainty in the parallax; with a perfect parallax, the uncertainty in 
the relative orbit produces a mass with an uncertainty of 0.9% (compared with the current value of 
2.5%), where the uncertainty in the relative orbit contributes equally to the total uncertainty in the 
mass. 

At that time this analysis can be done with considerable more precision, if our estimates of the 
temperature, and of the visual magnitude, or the absolute flux at any_ wavelength, of Sirius B, can be 
improved significantly over (he present values. At the moment, the uncertainty in parallax, 
temperature, and visual magnitude all contribute roughly equally to the uncertainty in the radius; in 
the future, the magnitude and temperature will dominate. This suggests that continued efforts to 
improve the accuracy of particularly the visual magnitude and of the temperature of the companion 
will be worthwhile. 

This work has been supported by NASA (grant NAG 8-564), by the NSF (grant AST-85-15747), by 
the A.J. Cannon Fund and by the Augustittds Fond in Denmark. 
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T efr\* 

5405 A 

1363 A 

30,000 K 

3.41e-4 

l'.98e*3 

35,000 K 

4.545e-4 

3.116e-3 

[H(35)-H(30)]/H(30) 

0.333 

0.574 

4T=1000 K 

AH=0.067 

4H=0.11'5 


AR=0.033 

^R=0.057 


Table 1: Showing the sensitivity of HfT^) to changes in T rff 

"■ '■ in the V band and the IUE spectral range.' “The sensitivity is least 
in the V band. From unpublished models by Shipman. 
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T _ 

eff 

HS 

Wesemael et al. 

HS log(g)— 8.5 

20,000 K 

1.70e-4 

1.69e-4 

(+0.005) 

30,000 K 

3.41e-4 

3.421e-4 

(-0.0082) 

50,000 K 

6.49e-4 

6.58e-4 

(+0.013) 


Table 2: Showing H , as calculated by Shipman and Wesemael et 

5405 A 

al. using completely different programs, for different and log(g). The 
first two coloumns are for log(g)=8 and the last is for the change in H when 
going from a log(g)=8 to a log(g)=8.5 model, relative to the log(g)—8 model. 
Note the very smaU dependencies of H on log(g) and choice of program. 
Modified from Table 1 of Wesemael et al. 1980. 
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Fig.l. Mass-Radius diagram for a 1.05 M 0 white dwarf, like Sirius B. The 

line marked C 12 is the Hamada-Salpeter zero temperature mass-radius 
relation and the two horizontal lines represent the radii calculated from the 
upper and lower limits on T gII - set by Holberg et al. Astrometric data and 
the magnitude was taken from Gatewood and Gatewood. If 1 sigma shifts 
are introduced in the magnitude and the parallax the radius-line given by 
2SGGQK moves down to the dashed line, showing the influence of errors on 
the input data. The velocity scale gives the equatorial rotational velocity in 
km/s for a homogenous, spherical body, modified from data given by 
Hamrick and Ostriker. The two circles have radii 1 and 2 sigma, 
respectively. 



Fig.2. Nomogram showing how Ta,f. radius and eauatorial rotational 
velocity relate for a 1.05 M 0 white dwarf like Sirius B. The calculations 
assume a spherical homogenous mass distribution and use data from 
Hartwick and Ostriker, Gatewood and Gatewood and unpublished spectra 
of Sirius B from Greenstein, Oke and Shipman. 
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ABSTRACT 

A 28 houc continuous EXOSAT observation of the 0.52 day period detached 

binary V471 Tauri has detected a strong, pulsed soft X-ray flux front the 

white dwarf conponent, with a 9.25 minute pulse period, an amplitude of 

^2055, and a double peaked pulse profile. A residual soft X-ray flux from 

the K dwarf companion is seen during white dwarf eclipse at orbital phase 

0.0. Pronounced dips in the soft X-ray light curve occur at orbital phases 

0.15, 0.18, and 0.85. The dips may be correlated with the wrtt t Lagrangian 

points of the binary orbit. Smaller dips at phases near the eclipse may be 

associated with cool loops in the K star corona at heights of 0.2 stellar 

radii. The X-ray flux from the white dwarf is consistent with thermal 

models for a white dwarf photosphere with T^-a, 3500 0 K, log g = 8. 0-8. 5, 

log N[He/H] < -4.5, and log N h - 18.65 +✓- 0.2. The K star has a luminosity 

30 - 1 

of 1.5 +/- 0.3 x 10 ergs s in the 0.03-2.5 keV band. 
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I. INTRODUCTION 

V47 1 Tauri (BD+16°516) is a well Known 0.52 day binary containing a hot 
DA white dwarf and a K2V detached companion. As a member of the Hyades 

o 

cluster (Young and Capps 1971) it has a presumed age of 5 x 10 years and a 
distance of 49 +/- 1 pc (Vandenberg and Bridges 1984). Young and Nelson 
(1972) derived a sini = 1.1 x lo" cm, ^sin^ i = 79°5 +/- 2 ? 0 , n K = 0.8 M 

9 

n^= 0.8 H 0 , R* = 6.5 x 10 cm. The white dwarf mass and radius are accurate 
to within 15%, and are consistent with the Hamada-Salpeter (1961) mass- 
radius relation for white dwarfs with log g = 8. 2-8. 5. Nodel atmosphere fits 
to IUE low dispersion spectra have determined T^ ** 35000 K and log g«/8.0 
(Guinan and Sion 1984). The absence of He II, C IV, and Si IV lines in the 
IUE spectra at a sensitivity of 0.5 J equivalent width constrain the 
photospheric helium abundance to be log NfHe/H] < -3.0 and constrain the 
carbon and silicon abundances to be less than 10 1 solar. (cf. Henry, 
Shipman, and Wesemael 1985). The K dwarf's synchronous rotation is one of 
the most rapid Known for a nondegenerate star in a detached system. ( Hall 

! 7 

Van Buren, Charles, and Nason (1980) identified V471 Tau as a HEAO A-2 
LED X-ray source. They considered the K star corona to be the liKely X-ray 
source, but did not rule out accretion of a K star wind onto the white 
dwarf. V471 Tau was detected as an Einstein Observatory IPC source by Young 
at al . ( 1983), who also considered the K star to be the probable X-ray 
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source/ although their observations did not cover white dwarf eclipse 
phases . 

In this Letter we report new results obtained iron a 28 hour continuous 
observation of V471 Tauri with the European Space Agency EXOSAT satellite t 
the detection of soft X-ray fluxes from both the white dwarf and the K 
dwarf/ the discovery of a 9.25 minute pulsation from the white dwarf/ and 
the discovery of orbital phase related soft X-ray dips . 

II. OBSERVATIONS 

The EXOSAT X-ray satellite was pointed at V471 Tauri continuously from UT 
1551 on 22 August/ 1985 to UT 1954 on 23 August, 1985, covering 2.24 binary 
orbits, including two white dwarf eclipses. During this time, 26.9 hours of 
data were obtained by the CMA at the focus of the LEI telescope (Korte et 
al . 1981), which is sensitive in the 0.03-2.5 keV bandpass. The thin lexan, 
aluminum, thick lexan, and boron filters were used. V471 Tau was detected 
in all but the boron filter. A weak 2-6 keV flux was detected by the medium 
energy (ME) detectors. He defer a report on this harder X-ray flux to a 
forthcoming paper. 

a) Soft X-rays from Both Stars - The LE light curve for the entire 
observation is shown in Figure 1 . The decline in the soft X-ray flux during 
the white dwarf eclipse phases ^centered at phases 1.0 and 2.0^, indicates 
that most of this flux comes from the white dwarf. The eclipse light curve 
is more clearly illustrated in Figure 2, which also shows a residual soft X- 
ray flux from the X star. The sharpness of the X-ray eclipse is consistent 
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with the 1 minute interval between optical contacts (Beavers, Oespez, and 
Pierce 1979), The count rates observed or estimated from each star in each 
of the four filters are given in Zable 1 . 

b) Soft X-ray Dips -* Pronounced dips in the soft X-ray flux at orbital 
phases 0.15, 0.18, and 0.85 can be seen in the LE light curve (figure 1). 
The phase 0.85 event is the strongest, with essentially conplete cutoff of 
the white dwarf flux at the depth of the event. It occurs relatively 
unchanged for three successive binary orbits. The phase 0.15 event occurs 
during the first orbit, but is gone 12.5 hours later. The phase 0.18 event 
also occurs during the first orbit, and is greatly diminished 12.5 hours 
later. There is also a broader, more shallow flux decrease between phases 
0.62 and 0.83 during the second orbit. This decrease does not occur during 
the next orbit. The -eclipse light curve (Figure 2) shows a lower observed 
flux level before the first eclipse ingress (top). During the next eclipse 
ingress (bottom), the flux has returned to the 'clean' level. The egress 
light curves, also shown in rigure 2, are characterized by sharper dips 
occurring between phases 0.04 and 0.06. 

c) White Dwarf Pulsations - A Fourier analysis of the LE flux has 
revealed a 9.25 minute periodic variation (Jensen 1985). A power spectrum 
for data obtained for the 'clean* phases 0.22 - 0.62 is shown in Figure 3. 
There is more power in the first harmonic at 3.604 mHz th.an in the 1.802 mHz 
period. The inset in Figure 3 shows the pulse profile for the subset of 
'clean' data obtained with the thin lexan filter. The pulse profile is 
double peaked, accounting for the strength of the first harmonic. The count 
rate at pulse minimum is *»80X of the count rate at pulse maximum. 
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XII. DISCUSSION 

a) The X-ray Fluxes - The soft X-ray fluxes from white dwarfs have 

generally been successfully fit by photospheric models (Kahn et al. 1984; 

Heise 1985; Petra, Shipman, and Canizares 1986). The white dwarf rates 

listed in Table 1 were fit by white dwarf model atmospheres for 

- tff = 30000 K, 32500 K, and 35000 K, and log g = 8.0 and 8.5 to test the 

possibility that the soft X-ray flux from the white dwarf in V471 Tau is 

photospheric in origin. Details concerning the model atmospheres we used 

can be found in Petre, Shipman, and Canizares (1986). Since the 

photospheric temperature, gravity, and distance to the source are already 

approximately Known, the observed EXOSAT count rates for the three filters 

provide a good constraint on the chemical composition of the photosphere and 

the interstellar column density. To model the composition of the 

photosphere, we varied the helium abundance as log H£He/Hl = -5.3 to -3.7 in 

steps of 0 . 1 dex. A pure hydrogen atmosphere model was also tested. He have 

assumed a negligible opacity due to metals heavier than helium. Simulations 

-3 

show that metal abundances less than 10 solar produce a negligible ..opacity. 
(Kahn et' al . 1984). IUE observations of V471 Tau have not constrained metal 
abundances in this range. Until more sensitive UV spectra are obtained, we 
cannot rule out the possibility that photospheric metal opacity has an 
effect on the emergent soft X-ray flux. 

The results from our fits are summarized in Table 2. Successful fits are 
found for T = 35000 K and log g = 8.0 and 8.5. A = 32500 K, 

log g = 8.5 model almost fits, suggesting that a finer grid of models would 
show a slightly higher temperature to be within acceptable bounds. The 
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agreement between tha EXOSAT data and white dwarf nodal atmospheres, f or 
tamparatures and gravities consistent with UV observations, provides 
evidence that the soft X-ray flux fron tha whita dwarf in V47 1 Tau 
originates in a ~ 35000 K photosphere Which nay contain trace alanants of 
heliun and/or heavier netals . The successful photospheric models give a 

-It -2 <-1 

0.03-2.5 keV white dwarf flux of 3 +/- 1 x 10 ergs cn s , corresponding 

to a luninosity of ip +/- | x 10* ergs s . Because the white dwarf flux is 

pulsed, it is probable that a fit of a honogenous photospheric nodel to the 
pulse averaged count rates does not correctly constrain the 'average* 
photospheric parameters. Additional fits of nora complex photospheric 
nodels are in preparation. 

The K star has a migrating RS CVn type wave (Young et al. 1983, and 

references therein), which implies spotted regions and a high density of 

magnetically confined coronal plasma. Its corona is presumably similar in 

nature to those of other late type stars observed by the Einstein and EXOSAT 

Observatories, and describable in terms of thermal emission with kT in the 

range 0.3-2 keV (a. g. Golub 1983). The observed ratio of rates in the thick 

and thin lexsrn filters is consistent with temperatures in this range for 

16 -7 

interstellar column densities of 1-7 x 10 cm , which are also consistent 
with the white dwarf models.' The implied 0.03-2.5 keV K star flux is 
~6 x 10 **ezgs cm 2 s 1 , corresponding to a luminosity of 1.8 x 10 30 ezgs s'* . 

b) The Dips - The lines of sight to the white dwarf during the dips are 
illustrated in Figure 4. Material located along these lines of sight was at 
least 2-3 stellar radii from the K star. The dips at phases 0.62 to 0.75 
must be due to material located closer to the white dwarf than to the K 
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star. Figura 4 shows that the dips at phasas 0.15, 0.18, and 0.85 occurred 
whan tha triangular - (outer) - Lagrangian points of the binary orbit are near 
tha lina of sight to tha white dwarf. These unique phases are 0.83 for the 
Lagrangian point on tha ingress side (L4) and 0.17 for the point on the 
egress side (L5) (Kopal 1978). Possibly ejected material can collect at 
these potential minima for at least several binary orbits. Alternatively, 
the material could be in very large coronal loops anchored to the K star. 
Another possibility is that absorbing material is located at greater, 
circumbinary distances, Sion and Bruhweiler (1986) have suggested the 
presence of cool circumbinary material in this system. 

^ Absorption dips in neutron star low mass binary systems have been 
attributed to structures related to accretion disks (e.g. White and Mason 
1986). There is no evidence in the V471 Tau UV data for an accretion disk 
(Guinan and Sion 1984). 

The sharp dips between phases 0.04 and 0.06 on the egress side of the 
eclipse are likely to be due to material bound to the. K star .. -The^event at 
phase 6.052 during tha first eclipse (top), if due to material in a cool K 
star loop, would require a loop height of ~0.2 stellar radii. The 
variability we observe in the soft X-ray light curve near the eclipse 
therefore suggests that dense, cool loops with heights ~ 0.2 stellar radii 
are present on the K star, and are variable in density and/or temperature on 
a time scale shorter than 12.5 hours. The transient C IV absorption line 
observed by IUE near the eclipse phase (Guinan and Sion 1984) is consistent 
with the presence of such structures. 
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c) The Whit* Dwarf Pulsations - He consider the most likely cause of the 
pulsations to he either 1) the rotation of a non-hoi»ogenous photosphere/ 
with the inhomogeneity produced by nass accretion onto magnetic polar 
regions of the white dwarf* or 2) nonradial g node oscillations excited by 
an. instability in the photosphere. 

i) Rotation? - The double peaked pulse profile is suggestive of a bipolar 

structure in the white dwarf photosphere* which could be readily explained 

if it is accreting nass fron the K star onto nagnetic polar regions. The 

nass accretion rate should deternine whether the poles are bright or dark in 

-ii,5 ~i 

X-rays. For example* a rate of 10 M^yr could produce 12. S eV polar hot 

,S 

spots with a filling factor of 10 . A blackbody flux fron such spots 

would be consistent with the EXOSAT data* and is an alternative to the 

photospheric models* but the required accretion rate is high for a detached 

binary in which the mass transfer is through capture of the K star wind. 

-n -I 

Accretion rates < 10 yr would not produce polar hot spots bright enough 

to modulate the soft X-ray flux* but may supply sufficient helium or metals 
to produce dark spots. Sion and Starrfield (198*0 have suggested that the 
photospheric metal lines observed in the hot DA2 white dwarf Feige 24 are 
supplied in this fashion from the detached M dwarf companion. The white 
dwarf magnetic field strength does not have to be as high as the fields in 
magnetic cataclysmic variables to be able to channel these lower accretion 
flows onto the poles . 

ii) Oscillations? - The instability of white dwarf photospheres to 
nonradial g mode oscillations* excited by the K -mechanism in partial 
ionization zones* is well established for DA white dwarfs sufficiently cool 
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(T<ff < 12000 K) to have unstable hydrogen partial ionization zones* the 
ZZ Ceti variables (Winget et al . 1982a). The discovery of DB pulsators wit) 
Tttf ** 30000 K (e.g. GD358* Winget et al. 1982b* and PG 1 654+1 60 , Winget et al . 
1984) confirms the existence of a helium instability strip near this 
temperature. Given the extreme sensitivity of the temperature of the heliur 
instability strip to the efficiency of convection (Winget et al . 1983)* it 
is possible that DA white dwarfs with T#f ^ 35000 K and a stratified H/He 
atmosphere could be unstable to nonradial g mode oscillations driven by a 
helium partial ionization zone beneath a thin hydrogen surface layer. Cox 
and Starrfield (1986) and Kawaler (1986) have recently argued that nonradial 
g node oscillations in hot DA white dwarfs could be driven by nuclear 
burning . 

Additional EXOSAT observations of V471 Tau in February 1986 confirm the 
existence of the 9.25 minute pulsations. The pronounced dips* which were 
prominent in the August 1985 observation* were not evident in the February 
1986 observations. A more detailed analysis and discussion of all of the 
EXOSAT observations of V471 Tau will appear in a forthcoming paper. 
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FIGURE 1 - Light curve of the 0.03-2.5 KeV flux from V471 Tauri for th 

entire 28 hour observation. The time bins are 554.85s, equivalent to on. 

white dwarf pulse period. The data from the thick lexan and aluminum filter: 

are normalized with respect to the thin lexan filter rate (cf Table 1), am 

combined with the thin lexan filter data to produce a light curve 

approximately equivalent to that which would be obtained from a continuous 

observation with the thin lexan filter. The horizontal dashed line indicates 

the mean thin lexan count rate between orbital phases 0.22 and 0.62. The 

vertical dashed lines mark the phases of quadrature. The phases at which the 
-friAoju W" 

-o u t a s; Lagrangian points (L4 and L5) are in our line of sight to the white 

dwarf are indicated by the arrows. The white dwarf is eclipsed at phases 

1 . 0 and 2.0. 

FIGURE 2 - Light curve of the 0.03-2.5 keV flux from V471 Tauri for the 
two eclipses of the white dwarf. The time bins Are 115.6s out of eclipse 
and 554.85s during eclipses. The top figure shows the first eclipse. The 
bottom figure shows the second eclipse. All of the data shown were obtained 
using the thin lexan filter except for phases 0.9-1.01 during ^ the ..first 
eclipse, which is thick lexan filter data. This data has been normalized to 
a level common with the thin lexan filter data. The data are also 

normalized with respect to a pulse profile template, so that the short time 

scale variability is not due to the pulse variability. 

FIGURE 3 - Power spectrum of the 0.03-2.5 keV flux for all data collected 
between orbital phases 0.22 and 0.62. The data from the thick lexan and 
aluminum filters are normalized with respect to the thin lexan filter rate, 
and combined with the thin lexan filter data to produce the input time 
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series . The power is normalized with respect to the Poisson noise level. The 
Hyquist frequency is 86.51 milliHertz (mHz). The frequency resolution is 
0.01056 mHz. The higher frequency power (unshown) is completely consistent 
with a Poisson noise process. The two horizontal dashed lines indicate the 
minimum power .levels for which features exceed Poisson noise expectations at 
confidences of 50% and 99%. 

FIGURE 4 - Schematic of the V471 Tauri binary system, indicating the 
orientations of the five Lagrangian points, and their Roche equipotentials , 
with respect to the white dwarf and K star components. A mass ratio q=1 is 
used. The K star clearly does not fill its Roche lobe. The curved arrows 
indicate the sense of rotation of the system. The dotted lines mark the 

orbital phase range for which dips are seen in the 0.03-2.5 keV light curve. 

4n<i 

The lines of sight to the white dwarf at the outer Lagrangian phases (0.167 
and 0.833) are indicated by dashed lines, as is the eclipse phase^ 0.0. The 
shaded regions indicate the phases at which dips are seen. The dip near the 
L4 point is particularly deep and persistent, and is therefore more heavily 


shaded. 
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Temperature, rtl<Uua and rotational velocity of Sirius B. 

by Peter ThejU and Harry L. Shipman. 

Physics Department, University of Delaware, 

Newark, Delaware 19716. 

Abstract. 

Based on the theoretical work, by Ostriker and Hartv/ick (1968), concerning the magnetic fields of 
rotating white dwarfs, and the recent determination of the effective temperature, by Holberg et al. 
(1S84), we have derived the interrelationships between radius, rotational velocity «n,d effective 
temperature for Sirius B. 

We find that i) the radii, 0.76 to 0.81 in hundredths of solar radii, predicted by T #R between 26000 
and 28000K and a mass of 1.05 M q places Sirius B somewhat above the Hamada-Salpeter mass-radius 
relation, by about 1 or 2 a. ii) the rotational velocity of Sirius B is predicted to lie between 0 and 
600 km/s (using 2 cr limits) and that only radii which place a white dwarf above the mass-radius line 
are allowed, iii) in order to further limit the rotational velocity and the radius of Sirius B better 
determinations of and magnitude are needed. If the satellite HIPPARCOS can determine the 
parallax of Sirius B with 1 milliarcsecond accuracy then the errors on T fff and the magnitude ' will be 
the limiting factors in the future. 

Keyword*: Sirius B • rotating white dwarfs - magnetic fields - Mass /Radius relation. 

L Introduction. 

The motivation for our work is twofold. We want to check if Sirius B falls on the mass-radius (M- 
R) relation using new data from Holberg et al. (1984). Only a few objects have so far been placed 
conclusively on that relation. We also wish to present a new method for determining rotational 
-velocities for cases- when the standard method, of modeling rotationaily broadened -lines, doesn’t work; - - 

Many determinations of Sirius B’s properties have been made over the years. Greenstein, Oke and 
Shipman (1971) analyzed the H q and H^ line profiles and arrived at T eff = 32000 +/- 1000K, along 

with gravity and radius: R/R 0 — 0.0078 +/- 0.0002. Shipman (1976) pointed out that a pure 

Hydrogen model with T #ff = 32000K would have the low absorption coefficient in the X-ray region 
necessary to explain the X-ray observations by Mewe et al. (1975). This value of T eff was obtained 
using unblanketed model atmospheres. If blanketed atmospheres had been used, this value of T^ 
would have been lowered by 2000K. Later determinations of T^ ff have confirmed the general 
conclusion that Sirius B is much hotter than Sirius A, but also lowered the T „ somewhat. Bohm- 
Vitense, Dettman and Kapranidis (1979) found T = 26000 +/- 1000K by fitting theoretical models to 
IUE spectra of Sirius B. Later Martin, Basri and Lampton (1982) used HEAOl A-2 data and fitted a 

28000K model finding support for a photospheric origin of the X- rays. Using the IUE data and UV 

data from Voyager 2 Holberg, Wesemael and Hubeny find support for T^ in the range 28000 • 28000K 
with a firm upper limit at 28000K. EXOSAT data will also provide a constraint on the T^, although 


1 


uncertanties in the EXOSAT calibration are still being investigated. 

One open question is whether, given all these determinations of T #ff , Sirius B fits on the M-R 
relation. One reason Sirius B might not lit the mass* radius relation is that rotation can increase its 
radius (Ostriker and Hartwiek, 1968). While one can analyze the rotationally broadened line profiles 
(Pilaehowski and Milkey, 1984) this method is only applicable to stellar spectra with very high 
resolution * and thus only useful for a few bright stars. Another possibility is that residual nuclear 
burning could increase the radius from its zero temperature value, as has been shown by Kocster and 
Schonberner (1986) for 0.6 M q stars. A similar calculation for the 1 M q case has not been done, but 
the effects of nuclear burning would be smaller, and possibly negligible, for a higher mass. 

Ostriker and Hartwiek established the relationship between radius, magnetic field and angular 
momentum for a white dwarf of mass 1.05M . We can use their work to get the relationships between 
radius, rotational velocity and temperature, provided that we have an upper limit to the magnetic field 
strength, and have values for orbital elements, the magnitude and the parallax. We have applied this 
method to Sirius B. 

Q. Method. 

From Keplers third law one can calculate the total mass of the Sirius system and use the accurately 
known t fractional mass of Sirius B to find its mass. The radius of Sirius B is found by using the 
calculated relationship between absolute flux and effective temperature (Shipman 1979, table 2) and the 
magnitude and parallax (Gatewood and Gatewood, 1978). In using <he Ostriker and Hartwiek (1968) 
results we assume that Sirius B rotates as a homogeneous sphere and that there is essentially no 
magnetic field on Sirius B. We base this last assumption on inspections of high-resolution spectra of 
the and lines (Greenstein, Oke and Shipman* 1971, unpublished spectra). 

As we used Ostriker and Hartwicks results we found that the v =0 line fell above the Hamada- 
Salpeter mass-radius relation. We assume this is because Ostriker and Hartwiek did not sufficiently 
model the envelope of white dwarfs and that if thev had done so the v =0 line would have fallen on 
the Hamada Salpeter M-R relation. Thus we put v^=0 on the M-R line for a C 12 body. We assumed 
the angular momentum of Sirius B could be modeled by using the moment of inertia of a homogeneous 
sphere. Using a more realistic mass distribution (/>(r)r 2 =const.) makes only a 20% difference in the 
equatorial rotational velocity v . 

In fag.1. we plot the results of our analysis in a M-R diagram. We have drawn lines parallel to the 
Bamada-Salpeter relation through each point determined by velocity and mass - we are thus assuming 
that the M-R lines for non-zero velocities will look like this in a small area close to the HS line. Fig.2 
is a more detailed nomogram where the relationships between T^, radius and rotational velocities can 
be seen for a 1.05M o , homogeneous, spherical white dwarf. 


2 


TIT- Uncertainties. 

The quantities that enter our calculations are of course in error to some extent - to see how 
important these errors are we both calculated the propagation of errors and introduced 1 a shifts in the 
parallax and the magnitude to see the effects. (See fig. 1.) The sixes of the relative errors on the 
quantities that determine the mass of Sirius B ate: 

W)- (*»>»* + + f'j/ + 

or 

(0.0261) 2 = (0.0246) 2 + (0.0079) 2 + (0.0030) 2 + (0.0022) 2 

where x,a,f_ and P are; parallax, angular semi-major axis of the orbit, the fractional mass of Sirius B 

D 

and the period of the orbit. The 3% error is mainly due to the error on the parallax, a contributes 
about one tenth as much. 

The error on the radius, determined from H(T <fJ ), r and magnitude m, is dominated by the error on 
the magnitude: 

(^..(Rb)) 2 - (tolO/S m ejmtf + (.„,(*))* + (. JH(TJ/2)f 

or 

(0.0192) 2 = (0.0166) 2 + (0.0082) 2 + (0.005) 2 

where we have assumed that the relative error on is 1% - a reasonable number, as will be 

discussed in section IV. The error on the magnitude contributes about four times as much as does the 
error on the parallax. We have used an absolute error on the magnitude of (0.03) +(0.02) where the 
last term takes into consideration the absolute calibration error. It is obvious that an improved 
determination of the magnitude is called for. Perhaps it is impossible, at the present epoch of the 
Sirius B orbit, and 'with standard techniques, to improve significantly on the value for the magnitude, 
but, as Smith and Terrile (1984) point out in their paper on the Pictoris disk the sort of equipment 
they used is very well suited for photometric work near bright sources, such as the Sirius system. 

How will the capabilities of the HIPPARCOS satellite change matters? If it can determine parallaxes 
to within 1 milliarcsecond then the err&r on Mg that z introduces will be the same as that from the 
semi-major axis and furthermore the uncertainty on the magnitude will overshadow that on the 
parallax by a factor of nearly 40 in its influence on the relative error of R Q . 

— r*™" TT— — " ■ ' — — — 
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IV. Discussion of the methods. 

We determine the radius R of Sirius B by using the fundamental relation between the stellar flux f, 
the monochromatic Eddington flux at the stellar surface H(T #ff ), and the stellar distance D, from 
Shipman 1979: 

(1) E = D (f/(4<rH(T "Vf* 

Because of the sensitivity of the results we discuss to relatively small uncertainties in the values of 
various input quantities, we include a considerable discussion of the above equation and the 
uncertainties therein, and ask for the readers’ patience with what is in places a pedagogical discussion. 
Earlier discussions of errors have been given by Shipman (1979) and, with respect to the accuracy of 
the H(T #ff ) relation, by Wesemael et aL (1980). * 

An advantage of this method, and of all surface brightness methods for determining a stellar radius, 
is that you do not need to know (or guess) how much energy the star is emitting at all wavelengths. 
The disadvantage is that you must be able to determine the stellar flux H accurately. But this 
problem can be circumvented or minimised by carefully choosing the wavelength or wavelength range 
over which you fit the above equation to the data. Provided that you apply the equation at a 
wavelength which is shortwaid of the Planck maximum, the dependence of H on is relatively weak 
- linear in the black-body approximation • which minimizes the uncertainty in R produced by 
uncertainties in T 

In this paper we apply it in the V band, since the V flux from Sirius B is known to within 0.03 
magnitudes. In contrast, Bohm-Vitense et al. applied the equation in the IUE spectral range, which 
one might think to be superior, being closer to the Planck maximum at these temperatures. However, 
model calculations indicate otherwise. H, _ increases by 33% when you increase T „ from 3O,000K 
to 35,000 K, while H increases by 57% (Shipman, unpublished, blanketed, non-LTE models^ the 

Wesemael et al. models are virtually identical^ as is shown in table 1). Thus, a 1,000 K uncertainty 
in T eff produces an uncertainty of 6.6% in H (and 3.3% in R) when you use this method in the V 
band, and an uncertainty of 11% in H (5.5% in R) when you use the method in the IUE spectral 
range. Thif consideration, by itself, suggests that the visual band is the one to work with. 

The ultraviolet spectral band is not as useful for application of equation (1), although it is essential if 
one is to determine T #ff accurately. The correction for scattered light in the IUE spectral range is 
considerable, being smallest at the shortest wavelengths, and is uncertain, since one has to use an 
existing spectroscopic instrument rather than a specially constructed area scanning photometer. In 
addition, the uncertainty in the overall IUE absolute calibration is usually quoted at 10% (see, e.g., 
R.C. Bohlin, et aL, 1980 and R.C. Bohlin, 1985). Individual leaders can make their own judgments 
of the accuracy of the IUE measurements of the Sirius B UV flux, but it is at least as uncertain as 
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the overall calibration uncertainty of 10%. Our examination of the figures in Bohm-Vitense et al. 
suggests that the uncertainty in measuring f generated by the need for correcting for scattered light is 
probably half of the total correction, or, at the shortest wavelengths, about 15%, producing an overall 
error in f of (10) + (15) or 18%, which corresponds to 9% in R arising from uncertainties in the flux 
measurements alone, if you apply equation (l) in the IUE spectral range. 

An additional consideration which must be applied in an investigation of this kind is how good are 
the models; how well can you trust the H(T gff ) relation which enters equation (1) ? This point has 
been discussed earlier by, <;.g., Auer and Shipman (1977), Shipman (1979b) , and most extensively by 
Wesemael et al. (1980). The ATLAS program, used for the H(T <ff ) relation used here (and in 
Shipman 1979), ;ind the Auer-Mihalas program as modified by Wesemael, used for the Wesemael et al. 
calculations, are completely independent. The uncertainties in H g405 for a given T^ are really quite 
small (see table 2.); at T eff = 20,000 K, comparing blanketed models with blanketed models* the 
relative change in calulated as (^ Shipmaa - H w «. ni «l/Hshi pm J “ +0-0059, at T # = 30,000 K it 

is <0.0032 and it is <0.0061 at T (ff — 50,000 K. Using unpublished, non- LTE models calculated at 
Delaware indicates a realtive change in H g405 A calculated as (^ nonLTE * ^lxe^^^LTE “ 0-011 at 
35,000 K. This scale is quite gravity independent; going from log(g) of 8 to 8.5 at 30,000 K decreases 
the Eddington flux, in these same relative terms, by -0.0082. The uncertainties in the H(T^) relation 
are thus far smaller than that of the magnitude and can safely be estimated at about 1% or less. 

V. Results. 

On Eg. 1 one sees that the temperature range that Holberg et al. arrived at corresponds to a range 
in radii of from 0.76 to 0.81 in 100R/R q . This is above the C mass-radius relation of Hamada and 
Salpeter, but if la shifts are introduced in magnitude, mass and parallax then Sirius B moves onto the 
M/R relation. Not until the errors on magnitude, parallax and of course T^ ff are reduced, as discussed 
above, can the issue of Sirius* B’s placement on or off the mass-iadius relation be settled. 

From the Holberg et al. temperatures and the Ostriker and Hartwick results we firfd the velocity to 
lie between 0 and 600 km/s - these are the velocity limits that a 2 a circle, centered at (M,R), touches. 
Furthermore, only positive rotational velocities are allowed, so only the part of the M-R plane above 
the C line is accessible. Any temperature and mass determinations that place Sirius B, or any other 
white dwarf, below the line can thus be ruled out. 

VL Discussion. 

Some interesting considerations, and some suggestions for current work, emerge from a consideration 
of how the astrometric data on Sirius are likely to improve in the future,, The European satellite 
HEPPARCOS ean improve the parallax of Sirius significantly; the current best parallax (Gatewood and 
Gatewood, 1978) has an uncertainty of Z milliarcseconds, and the advertised B1PPARCOS capability 
should improve this by a factor of 3. Such improvement can, in principle allow us to determine the 
radius with much greater precision than we can now. In addition, at the present time the uncertainty 
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in the mass is dominated by the uncertainty in the parallax; with a perfect parallax, the uncertainty in 
the relative orbit produces a mass with an uncertainty of 0.9% (compared with the current value of 
2.5%), where the uncertainty in the relative orbit contributes equally to the total uncertainty in the 
mass. 

At that time this analysis can be done with considerable more precision, if our estimates of the 
temperature, and of the visual magnitude, or the absolute flux at any_ wavelength, of Sirius B, can be 
improved significantly over the present values. At the moment, the uncertainty in parallax, 
temperature, and visual magnitude all contribute roughly equally to the uncertainty in the radius; in 
the future, the magnitude and temperature will dominate. This suggests that continued efforts to 
improve the accuracy of particularly the visual magnitude and of the temperature of the companion 
will be worthwhile. 

This work has been supported by NASA (grant NAG 8-564), by the NSF (grant AST-85-15747), by 
the A.J. Cannon Fund and by the Augustinus Fond in Denmark. 
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T « ff \ x 

5405 A 

1363 A 

30,000 K 

3.41e-4 

L98e-3 

35,000 K 

[H(35)-H(30)]/H(30) 

4.545e-4 

3.116e-3 

0.333 

0.574 

4T=1000 K 

4H=0.067 

^1H=0.115 


4R=0.033 

AR=0.057 

...Table JLi . 

Showing the sensitivity of H(T eff ) to changes Jn 

T 


in the V band and the IUE spectral rang*. The sensitivity is least 
in the V band. From unpublished models by Shipman. 


T 

eff 

HS 

Wesemael et al. 

HS log(g)=8.5 

20,000 K 

1.70e*4 

l,69e-4 

(+0.005) 

30,000 K 

3.41e-4 

3.421e-4 

(-0.0082) 

50,000 K 

6.49e>4 

6.58e-4 

(+0.013) 


Table 2: Showing H g405 A , as calculated by Shipman and Wesemael et 

ah using completely different programs, for different T eff and log(g). The 
first two coloumns are for log(g)=8 and the last is for the change in H when 
— - going from -a log(g)=8 to a log(g)=8.5 model, -relative to the log(g)=8 model. 

Note the very small dependencies of H on log(g) and choice of program. 
Modified from Table 1 of Wesemael et aL 1980. 
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Fig.l. Mass-Radius diagram for a 1.05 M 0 white dwarf. like Sirius B. The 

line marie d C 12 is the Hamada-Saipeter zero temperature mass-radius 
relation and the two horizontal lines represent the radii calculated from the 
upper and lower limits on T el -f set by Holberg et al. Astrometric data and 
the magnitude was taken from Gatewood and Gatewood. If 1 sigma shifts 
are introduced in the magnitude and the parallax the radius-line given by 
2SG00S moves down to the dashed line, showing the influence of errors on 
the input data. The velocity scale gives the equatorial rotational velocity in 
km/s for a homogenous, spherical body, modified from data given by 
Hartwick and Ostriker. The two circles have radii 1 and 2 sigma, 
respectively. 





Fig.2. Nomogram showing how T gl f, radius and equatorial rotational 
velocity relate for a 1.05 M 0 white dwarf like Sirius B. The calculations 
assume a spherical homogenous mass distribution and use data from 
Hartwick and Ostriker, Gatewood and Gatewood and unpublished spectra 
of Sirius B from Greenstein, Oke and Shipman. 
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